A complete CO 2-1 map of M51 with HERA 
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Summary. The nearby, almost face-on, and interacting galaxy M51 offers an ex- 
cellent opportunity to study the distribution of molecular gas and the mechanisms 
governing the star formation rate. We have created a complete map (Fig. 1) of M51 
in 12 CO 2-1 at a resolution of 11" corresponding to 450 kpc using HERA at the 
IRAM-30m telescope. In [1] we have combined these data with maps of HI and the 
radio-continuum to study the star formation efficiency, the local Schmidt law, and 
Toomre stability of the disk in radial averages out to radii of 12 kpc. Here, we also 
discuss the distribution of giant molecular associations and its mass spectrum, in 
comparison with similar studies in the literature. 



1 CO, HI, and the radio continuum 

M51 is an interacting, grand-design spiral galaxy at a distance of only 8.4 Mpc 
seen nearly face-on. The emission detected with the 30m telescope (Fig. 1) 
traces the well known two-armed spiral pattern out to the companion galaxy 
NGC5195, which shows up brightly in the north-east at ~ 10.5 kpc radial 
distance, and out to the south-western tip of the second arm at the opposite 
side of M51. The outer parts of the two arms in the west and in the east appear 
more fragmented than the inner parts. The western arm especially is almost 
unresolved. Inter-arm emission is detected above the 3cr level out to radii of 
about 6 kpc. Several spoke-like structures connect the spiral arms radially. 
The large-scale distribution of the 21 cm line of atomic hydrogen in M51 was 
analyzed by [2] using the VLA. The HI emission at 13" resolution (Fig. 2) is 
weak in the inner region while the outer CO arms are clearly delineated in 
HI. These data have been combined with a new large scale 20cm map of [3] 
which provides an extinction-free estimate of the star formation rate. 
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Fig. 1. Map of 12 CO 2-1 integrated intensities in Kkms" 1 showing M51, i.e. 
NGC 5194 and its companion galaxy NGC 5195 in the north-east. Crosses mark 
the center positions of the identified 155 giant molecular associations (GMAs). 

We derive a global star formation rate of 2.56 M Q yr -1 . The total gas 
surface density J2 ga .s = 1-36(X)h 2 +Shi) drops by a factor of ~ 20 from 
70 M Q pc~ 2 at the center to 3 M Q pc -2 in the outskirts at radii of 12 kpc. The 
ratio of HI over H2 surface densities, X)hi/Sh~2> increases from ~ 0.1 near 
the center to ~ 20 in the outskirts without following a simple power-law. The 
star formation rate per unit area drops from ~ 400 M Q pc -2 Gyr -1 in the 
starburst center to ~ 2M Q pc -2 Gyr -1 in the outskirts. X] gas an d Ssfr are 
well characterized by a local Schmidt law X)sfr x Sgas with a power-law 
index of n = 1.4 ±0.6 (Fig. 3). 

The critical gas velocity dispersions needed to stabilize the gas against 
gravitational collapse in the differentially rotating disk of M51 using the 
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Fig. 2. VLA map of integrated HI intensities [Jy/beam] at 13" resolution [2] in 
colors. Contours show integrated 12 CO 2-1 intensities (cf. Fig. 1). 

Toomre criterion, vary with radius between 1.7 and 6.8 kms" 1 . Observed radi- 
ally averaged dispersions derived from the CO data vary between 28kms~ 1 in 
the center and ~8kms _1 at radii of 7 to 9kpc. They thus exceed the critical 
dispersions by factors Q gas of 1 to 5. Taking into account, in addition, the 
gravitational potential of stars, the disk may be critically stable. 

2 Distribution of molecular gas 

Molecular clouds are the sites of all star formation in the Milky Way and 
also in external galaxies. It is therefore of great interest to study the global 
distribution of star formation in entire galaxies by studying the distribution 
of the molecular gas. Nearby face-on galaxies like M51 offer the possibility 
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Fig. 3. Radially averaged star formation rate per unit area, X/sfr' versus surface 
densities of the total gas X/ ga s' an< ^ °^ ^ 2 an< ^ ^1 ori fy- The solid green line is the 
local Schmidt-law found in M51. The dashed-dotted black line is the global Schmidt- 
law found by [4]. Drawn blue lines represent lines of constant gas depletion time or 
star formation efficiency. 

to study the distribution of molecular gas over their entire surface at high 
spatial resolution and without the distance ambiguities as encountered for 
the Milky Way. In M51, the total mass of molecular material derived from 
the integrated 12 CO 2-1 intensities is 21O 9 M0. The 3er limit with resolu- 
tions of 11" and 5kms _1 corresponds to a mass of 1.710 5 M Q . The spatial 
resolution of 450 pc resolves structures larger than typical GMCs which we 
label giant molecular associations. These may be bound clusters of GMCs as 
suggested by [5] or random superpositions of GMCs [6] . In a first attempt to 
study the properties of these GMAs, we have decomposed the CO 2-1 data 
set into Gaussian shaped clouds using the gaussclumps algorithm [7]. This 
algorithm has been developed to study the statistical properties of the struc- 
ture of molecular clouds as seen in spectral line emission cubes and has been 
applied to many Galactic clouds, e.g. [8, 15]. It iteratively fits a local Gaussian 
to the global maximum of the data and subtracts it. The algorithm fits the 
intrinsic, deconvolved widths to the observed data taking into account the 
given angular and spectral resolutions. It outputs the clump center positions 
and widths in the two spatial and in the velocity coordinates. It also outputs 
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the clump peak temperatures, as also their orientations. In the following, we 
will use the term clouds or GMAs rather than clumps. 
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Fig. 4. The mass spectrum of GMAs identified in M51 shows the number of sources 
per logarithmic mass interval. Errorbars represent \j AN statistical errors. The best 
fitting power law of the form dN/dM oc M~ a has a slope of a — 2.0 ± 0.2. 



Table 1. Mass spectra of GMCs and GMAs in external galaxies. The number 
of identified clouds is given in column 4. M m j n and M max are the minimum and 
maximum cloud mass detected. Mt urn is the mass at the turnover of the spectrum. 



Source Distance Resol. No. a M m i n M tur n Af max Lit. 

[pc 2 ] [M ] [M ] [M ] 



NGC 4038/39 


19 Mpc 


310 x 480 


100 


1.4 


± 


0.1 


2 


10 6 


5 


10 6 


9 


10 8 


[9] 


M51 


8.4 Mpc 


450 


155 


2.0 


± 


0.2 


5 


10 s 


3 


10 6 


1 


10* 


[10] 


M33 


850 kpc 


50 


148 


2.6 


± 


0.3 


3 


10 4 


2 


10 6 


7 


10 6 


[11] 


M31 


780 kpc 


90 


389 


1.6 


± 


0.2 


2 


10 4 


1 


10 5 


5 


10 5 


[12] 


LMC 


54 kpc 


41 


168 


1.9 


± 


0.1 


4 


10 4 


8 


10 4 


3 


10 6 


[13] 



In M51, the algorithm decomposes 78% of the total mass into 155 clouds. 
Their deconvolved, i.e. intrinsic, sizes and FWHMs are larger than 20% of the 
resolutions. The 16 most massive GMAs with masses between 2 10 7 M Q and 
1 10 8 M follow the two inner logarithmic spiral arms like beads on a string 
(Fig. 1). These masses exceed the masses derived by [5] using the three-element 
Caltech Millimeter Array by a factor of ~ 2, presumably due to the slightly 
larger beam size of the present study and missing short spacing information. 
Such very massive molecular cloud complexes have also been found in other 
galaxies. See the compilation by [9] (their Table 3). In our study of M51, 
GMAs are also identified along the two outer arms extending towards the 
companions galaxy in the north and to the south. At least five spoke like 
features connecting spiral arms are decomposed into clouds. Only few clouds 
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are seen in the inter arm medium. The least massive clouds have a mass of 
510 5 M Q . 

Figure 4 shows the mass spectrum of all 155 clouds. It follows a linear 
slope of a = 2.0 ± 0.2 above the turnover mass of 5 10 6 M Q . The turnover is 
presumably not intrinsic to the GMA distribution but caused by the detection 
limits, similar as in Galactic clouds [8]. A slope of 2.0 agrees well, within the 
error, with more recent large-scale mapping studies of Galactic clouds by e.g. 
[14, 8, 15, 16]. These studies find a rather constant slope of a = 1.6 — 1.9 for 
cloud and clump masses between more than 10 4 M Q and 1O _4 M . Todate, 
only few studies exist which derive mass spectra of the entire cloud population 
of external galaxies. In Table 1, we compare the mass spectrum of M51, with 
those found in NGC 4038/39, M33, M31, and the LMC. The slope which we 
find in M51 is similar to the slope found in the Andromeda Galaxy using the 
same algorithm, though the study of M31 covers a lower mass range due to 
its proximity. The slope is also similar to the slope found in the nearby dwarf 
galaxy LMC. The Antennae exhibit a rather flat mass distribution with very 
massive GMAs while M33 shows a very steep spectrum with a = 2.6 similar 
to the initial mass function. 
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